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C Abstract 

' We construct a family of chaotic dynamical systems with explicit broad dis- 

tributions, which always violate the central limit theorem. In particular, we 
\ show that the superposition of many statistically independent, identically dis- 

\ tributed random variables obeying such chaotic process converge in density 

' to Levy's stable laws in a full range of the index parameters. The theory 

related to the connection between deterministic chaos and non-Gaussian dis- 
. tributions gives us a systematic view of the purely mechanical generation of 

Q\ ', Levy's stable laws. 
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The central limit theorem (CLT) breaks down when any stationary stochastic process has 
infinite variance. Thus, it is an important question whether this mathematically pathologi- 
cal situation is physically relevant or not. One noteworthy point arising from recent studies 
is that such examples of Levy's stable laws, which are the most famous class of distributions 
violating the CLT, can be seen in many different fields, such as astronomy, physics, biology, 
economics and communication engineering, under broad conditions Thus, it is our 

primary interest to determine why such Levy's stable laws are widely observed, and fur- 
thermore, to elucidate the mechanism of generating Levy's stable laws. In the 1980's, there 
were several studies which clarified the relation between intermittent periodic mapping and 
anomalous diffusion with Levy's law-like broad distributions . Random- walk models ^ 
and combinations of several random number generators are also utilized to generate 
Levy's stable laws. However, such analyses include approximations or the non-deterministic 
nature in the models themselves, or else their generation methods are only applicable to a 
special class of Levy's stable laws. The purpose of the present paper is to present a sys- 
tematic method for exact and purely mechanical generation of stable laws with arbitrary 
indices only using concrete chaotic dynamical systems. Let us consider an one-dimensional 
dynamical system 

^n+l = ^(X„-^) = /(X„) (1) 

on the infinite support (—00, +00). Note that this mapping f{X) can be seen as the doubling 
formula of — cot(^^) as — cot(26') = /[— cot(6')]. Thus, the system has the exact solution 
Xn = — cot(|2"6'o). Using a diffeomorphism x = 4>^^{9) = — ^^^^hq) of ^ G [0,2] into 

] — 00, +00 [, we derive the piecewise-linear map g^'^\0) = o / o (j)^^(0) as 

g(^\9) = 2e-2, 9e[l,2]. 

Because the map (0) has the mixing property (thus, is clearly ergodic) and preserves the 
Lebesgue measure ^d9 of [0, 2], the map / preserves the measure 

This is an explanation of the mechanical origin of the Cauchy distribution (|^). Note that the 
Cauchy distribution is a simple case of Levy stable distributions with the characteristic a = 
1. The measure @ is absolutely continuous with respect to the Lebesgue measure, which 
implies that the Kolmogorov- Sinai entropy h{fi) is equivalent to the Lyapunov exponent of 
ln2 from the Pesin identity; and the measure is a physical one in the sense that, for almost all 
initial conditions Xq, the time averages lim^^ooll/"^) J27=o ~ reproduce the invariant 
measure fi{dx) Our next step is to generalize exactly solvable chaos ([^) to capture the 
full domain of Levy's stable laws. Now, let us consider the mapping 

X„+i = |^(|X„r - l/l^nl")!* ■ Sgn[(X„ - 1/Xn)] = UXn), (4) 

where sgn(x) = 1 for x > and sgn(x) = —1 for x < 0. We prove here that this chaotic 
dynamics (P also has the mixing property similar to mapping (||), as well as the exact 
invariant density function given by 
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Pa{x) = —-, — . .„ . — — l^l for \x\ — >• oo. (5) 



Note that the chaotic system is a special case of with a = 1. We remark here that this 
system can also be seen as a doubling formula 5(26*) = /q,[s(6')], where s{6) = ~-j^^^^^tt^- 
Using the relations 

s{29) = U[sm for 9e[0,l) 
s{2e-2) = f4s{d)] for ^€[1,2] 

and defining the diffeomorphism 

^ ^ |tan(f^)|V« 

of 6* G [0, 2] into ] — oo, +oo[, we obtain the piecewise-linear map g'^'^\0) = (pa ° fa ° ^Pa^i^) 
as 

^ ^^^-\20-2, ^G[l,2] 
with the invariant measure ^d6 of [0,2]. Thus, map preserves 



/.(rfx) = p.(x)dx = -^dx = (9) 

Therefore, the class of dynamical systems (^) with the parameter a also has a mixing 
property (thus, is ergodic) with the Lyapunov exponent ln2. 

More generally, from the family of Chebyshev maps Yn+i = f{Yn) defined by the addition 
formulas of the form sin^(p^) = /[sin^(^)], where p = 2,3,---, with the unique density 
a{y) = — , ^ of the logistic map = 4F„(1 — y„) (which corresponds to the case p = 2) 



||TT| , [T2|| , we may construct infinitely many chaotic dynamical systems = /a(X„) with 

the unique density function Pa{x) given by Eq.(^ [|13[. For example, an explicit mapping 
with the Lyapunov exponent ln3 is given by 

V f (V \ \ - 3) ,i X„(|X„p°-3) , . 

X.^i = /.(X„) =1 (3|^^p„ _ 1) h -^^g-l (3|X„|--1) 

which has the density (^) from the triplication formula of s{9). 

In this case, the topological conjugacy relation g^^\9) = 0q o o (j)~^{9) yields the 
piecewise-linear map: 



e [0, 3, 

36-4, ee[l2]. 




g^'\9) = { 39-2, 9e[l,l) (11) 



In general, the same kind of topological conjugacy relation g^P\9) = 0„ o /„ o (p-\9) holds 
for a p-to-one piece- wise linear mapping g^P\9). Let us consider slightly modified dynamical 
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systems Xn+i = fa,s{Xn) = Ifai^^n) with a change of variable h{x) = Sx for a constant 
S > 0. Thus, this modified dynamics, 

faAX) = \li\Xr - l/\S'xn\^ . sgn[X - ^], (12) 

has an invariant measure Pa,s{x)dx = 6pa{Sx)dx = J^'|JjL|^|2a) with a shghtly modified 
power-law tail as 

a 1 

^(x) ~ — — — - for X ^ ±oo. (13) 

We will show that this power-law tail of density is sufficient for generating arbitrary symmet- 
ric stable laws. The canonical representation of stable laws obtained by Levy and Khintchine 
MM is: 



1 

P{x;a,P) = — / exp{izx)-ip{z)dz, (14) 
zvr J-oo 

where the characteristic function ip{z) is given by 

il^iz) = exp{— 272; — r7|;z|"[l -|- iPsgn{z)u!{z, a)]}, (15) 

a, P, 7 and i] being real constants satisfying 0<a<2, — 1</3<1,7>0 and 

07(2, a) = tan(7ra/2) for 07^!, 
(^^(z, a) = (2/7r) log 1^1 for a = 1. 

According to the generalized central limit theorem (GCLT), it is known WE\ that if the 
density function of a stochastic process has a long tail, 

p(a;) ~ c_ for x ^ — 00, p{x) ^ c^\x\^^^^°'^ ior x -|-oo, (17) 

then the superposition Sn = iJ2iLiX{i) — A]sf)/B]\f of independent, identically distributed 
random variables X{1), ■ ■ ■ , X{N) with the density p(x) converges in density to a Levy's 
stable law P{x;a,j3) with 

(3= {c+-c^)/{c+ + c_), 

AN = O,B^ = N^/'^,r^ = ^,,0±±^, for < « < 1, (18) 
An = N<x>,B^ = iW^ ^ = ^ , for l<a<2. 

In the case of the symmetric long-tail (|^), c+ = c_ = = and r) is determined by 



r] = ....,lc^■n,^^ for < a < 1 

^ ~ a5«sin{2^)r(a-l) 



ff^vFT — TT for 1 < a < 2. ^ ' 



Thus, according to the GCLT, the superposition of statistically independent, identically 
distributed random variables generated by unique chaotic systems (|) is guaranteed to 
converge in distribution to an arbitrary symmetric Levy's stable law P(x; = 0). Figure 
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1 shows that the convergence in distribution to a Levy's stable distribution with parameters 
a = 1.5 and /3 = is clearly seen for = 10000, as predicted by the GCLT. What 
about more general stable distributions including asymmetric stable laws? In the next part, 
we will also provide chaotic dynamical systems with explicit broad distributions, whose 
superpositions converge in distribution to arbitrary asymmetric Levy's stable laws. Let us 
consider a family of dynamical systems Xn+i = fa,Si,S2{^n), where 



fo. X>± 



(20) 



We can show that the dynamical system Xj+i = fa,5i,S2{^i) has an asymmetric invariant 
measure, fi{dx) = p{x;a,6i,62)dx, where Pa,5i,S2{.x) is given by 



Pa,Si ,62 (•^) 



a5fx 



a^a—l 



for X > 0, Pa,Si,S2i^) 



7r(l + 5i2"a;2") 
Because the power-law tail is asymmetric as 



a6^ 
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2a ~) 



for X < 0. (21) 
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(22) 



for 61 62, the GCLT guarantees that the limiting distribution would be a Levy's canonical 

^ 7^ 0. Thus, we can generate arbitrary 



(5?-<5? 



form P{x; a, (3) with the skewness parameter (3 

Levy's stable laws P{x] a, /3) for < a < 2 and — 1 < /3 < 1 [|l^ using the chaotic mappings 
/a,5i,52(^) with proper parameters a, 5i and 52- Convergence to an asymmetric Levy's stable 
distribution with indices a = 1.5 and [3 = — is clearly seen for = 10000 in Fig.2, 
as predicted by the GCLT. To show the exactness of the asymmetric density (^), we must 
check that the invariant measure Pcx,&^,&2i.^)'^^ satisfies the probabihty preservation relation 
( Perron- Frobenius Equation) ||18||: 



Pa,5uS2{y) 



~ fa, Si, 62 



X 



dx 
dy 



We note that fa,sA^) = faA^) and pa^M^) = Pa,Sii^) for a; > and Pa^M^) 
for X < 0, which also have the Perron-Frobenius equations 



(23) 



Pa, (52 ('^) 



paAiy) 



E 

=/«i(j/) 



dx 



for 



1,2. 



(24) 



Here, we define two pre-images Xa = fa}i,52^y) < and Xb = fa}i,52^y) > for y > 0. 
In the case fa^sA^) = y > 0, also define two pre-images x'^ = faSiiy) < and x'^ = 
fa5i{y){= ^b) > for ?/ = fa,5i{x) > 0. It is easy to check that 52Xa = 5ix'^. From the 



Perron-Frobenius equations (|23D and (|J), we have the relation 
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\X—Xa I (^3^ \X—X^ 



However, we can check the relation ( pSf ) under the condition (52Xa = Six'^. In the same 
manner, we can also show that paMM^v) satisfies the Perron- Frobenius equation (|23|) for 
y<Q. 

There is also an interesting dualistic structure in these types of chaotic dynamical sys- 
tems. Let us consider dynamical systems X*^^ = /* &2^-^n) defined as 



/a,5i,(52(^*) 



KkP^)^ for A-<-i 

Because the normalized and symmetric case, = /^,5j=i,52=i(-^*), can also be seen 

as the doubling formula of — sgn[tan(|6')]| tan(|^)|^/", and we can derive the piecewise- 
linear map g^'^\9) = (p*a ° fa ° (D by the diffeomorphism x* = (f>*a~^{9) = 

— sgn[tan(|6')]| tan(|6')|-^/", map f*{X*) also has the same invariant measure Pq,(x) given 
by Eq. (H). Thus, in the same way as that used in obtaining Pa,Si,S2{^) for fa,5i,S2{^)y 
we can show that map (|26|) has the invariant measure Pa,i/Si,i/S2{^)d^- Furthermore, the 
dynamical reciprocal relation 

f:,i/5,,i/S2iX*)faA,sM) = '^ for XX* = 1 (27) 

holds. This duahstic structure of the dynamical systems fa,Si,S2{^) and fa,i/Si,i/S2(-^*) 
originates from the relation (p~^{9) ■ (p^^{0) = 1. 

In summary, we have found that Levy's stable laws can be directly generated by the 
superposition of many independent, identically distributed dynamical variables obeying cer- 
tain unique chaotic processes. Owing to the ubiquitous character of Levy's stable laws, our 
methods related to ergodic transformations with long-tail densities have a broad range of 
applications to many different physical problems where stable laws play an essential role. 
This work is supported by RIKEN-SPRF and RIKEN-PSR grants. 
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FIGURES 

FIG. 1. (color) Densities of the superposition Sn = (Z]ili-^(0 ~ An)/Bn of dynamical 
variables X{i) generated by chaotic systems Xj^i{i) = /^^j^ 5[Xj(i)], with N different initial 
conditions -'^o(0|j=i, -,iv are plotted for N — 1,10,100,1000 and 10000. In this case, the 
limit density converges to the symmetric Levy's stable law with the indices a — 1.5 and 
= 0. 

FIG. 2. (color) Densities of the superposition Sj^ — (Eili^(0 - An)/Bn of dynamical 
variables X{i) generated by chaotic systems Xj^i{i) = fali.5Si=iS2=o.5[-^ji'^)]^ with N dif- 
ferent initial conditions Xo(i)\i=i,...^N are plotted for N = 1,10,100,1000 and 10000. In 
this case, the limit density converges to the asymmetric Levy's stable law with the indices 
a = 1.5 and f3 = -9^. 
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alpha=1 .5,beta=0[Fig. 1(Ken Umeno)] 
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alpha=1.5,beta=-frac{9-4sqrt{2}}{7} [Fig.2(Ken Umeno)] 



N=1 ^ 
N=10 -+-- 
N=100 




-10 -5 5 10 

X 



